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Abstract 

Spectroscopic observables such as electromagnetic transitions strengths can be related to the 
properties of the intrinsic mean-field wave function when the latter are strongly deformed, but 
the standard rotational formulas break down when the deformation decreases. Nevertheless there 
is a well-defined, non-zero, spherical limit that can be evaluated in terms of overlaps of mean- 
field intrinsic deformed wave functions. We examine the transition between the spherical limit 
and strongly deformed one for a range of nuclei comparing the two limiting formulas with exact 
projection results. We find a simple criterion for the validity of the rotational formula depending 
on (AJ" 2 ), the mean square fluctuation in the angular momentum of the intrinsic state. We also 
propose an interpolation formula which describes the transition strengths over the entire range of 
deformations, reducing to the two simple expressions in the appropriate limits. 
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I. INTRODUCTION 



In mean-field theories, electromagnetic transition rates are often evaluated using the rota- 
tional formula jl | to relate them to the multipole moments of the mean-field wave functions. 
The formula is justified by factorizing the wave function as a product of a wave function 
for the orientation angles times an intrinsic wave function and assuming that the matrix 
elements between intrinsic states at different orientations vanish. From a more microscopic 
point of view, the formula can be obtained as the strong deformation limit of the transition 
probability computed with angular momentum projected wave functions There are 



several studies in the 
deformed nuclei 
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iterature investigating the validity of the rotational formula in well 
ij]. However, as far as we know there has never been a systematic 
study of the validity and eventual breakdown of the rotational formula as the wave function 
approaches the spherical limit. A motivation for this study is the wide-spread use of this 
formula even outside of its domain of validity. For example, the increasing popularity of the 
Bohr Hamiltonian (| as a tool to handle low energy vibrational and rotational properties in 
a mean-field framework calls for a careful analysis of the limitations of the rotational formula 
for B(E2) transition strengths[7]. Often near-spherical configurations have a non-negligible 
amplitude in the wave functions and their contribution to the transition strengths needs 
to be handled with care. The purpose of this paper is to establish criteria for the use of 
rotational formulas, as well as to find useful approximations simpler than the full angular 
momentum projection to deal with moderate and soft deformations. 

This paper is organized as follows. Sect. II below discusses the representation of the wave 
function at small deformations. Our main result, derived in Sect. Ill, is an an expression for 
the transition strengths valid for small deformations Eq. (fTTj) below. This expression gives a 
non-zero value in the limit of vanishing deformation, in contrast with the rotational formula, 
Eq. ([3]) below. In Sect. IV we examine the validity of the formulas by comparing with full 
projections from the intrinsic states, taking a number of representative examples including 
quadrupole and octupole transitions. The dividing line separating the small and large de- 
formation limits is seen to be closely connected to the the angular momentum content of 
the intrinsic wave function. This gives a simple criterion to identify the regions of validity 
of the rotational formula. We also find that the B(E2) values can be simply parameter- 
ized as a function of the the quadrupole deformation parameter, Eq. (TT8j) below. Other 
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transition strengths like the B(E3, 3 _ — > + ) will be discussed and we will see that similar 
considerations apply to them as well. 

To set the notation, the multipole operators are defined as js] 



and the corresponding electric operators as 

Q%i = e^^Qx,. (2) 
The rotational formula for an axially symmetric intrinsic state is given by 

5( J EJ;J^O) R oT = ^|(0|g e Jo l0)r- (3) 

II. MEAN FIELD WAVE FUNCTIONS NEAR SPHERICITY 

The first step is the characterization of the intrinsic wave functions near sphericity. We 
will focus on quadrupole deformation because the generalization to other multipolarities 
is straightforward. We assume that the intrinsic wave functions are of the Hartree-Fock- 
Bogoliubov (HFB) mean-field type. The wave function \4>{q)) is labeled by the components 
of the quadrupole moment q 2fl = (4>\Q2n\4>) (/•* = ~ 2, ...,2). The wave function can be 
expressed in terms of a suitable spherical reference state |0(O)) by means of the generalized 
Thouless theorem 

10(g)) =jV ? exp(^(g)i0(O)>. 

Here Z{q) is a sum of 2-quasiparticle creation operators and Af q is a normalization constant. 
Given the Bogoliubov amplitudes U(q), V(q) and U(0), V(0) defining \<f){q)) and |0(O)) (see 
j^l for notation) the explicit form of Z(q) can be obtained p, App. E.3]. However, we only 
need to assume for the formal development below that Z can be expanded as a power series 
in q, 

z{q) = $> M (-ir4_, + 1 ^ Wv (-ir+^_^, + • • • 

The phases are introduced for consistency with the following properties of the defor- 
mation parameters q 2 ^ = (Q 2fJ ) = (^2^)* = (-1)^(^2-^) = (-l)^-^- It also im- 
plies that Z ( 2 l l + = (-iyZ 2 ,-» and zf^, = (-1)^' Z' 2 The tensor character 
of the multipole operators implies that the deformation parameters of the rotated wave 



function |0(<jO) — R(^)\4 > ( ( l2fj)) also behave as the components of a spherical tensor 
<?2^ = Yl/jj,' ^M^)5V- To be consistent with this property, the operator Z 2 ^ must transform 
under rotations as 

The corresponding transformation properties of the operators Z' 2 , are given by 

vv' 

This property makes it possible to decompose the operator as the direct sum of spherical 
tensors 

JM 

In the present example the range of the spherical tensors Z' JM is J = 0, ... ,4. Using the 
same kind of arguments it is easy to show that the Zand Z' operators must be even under 
parity. The generalization to an arbitrary multipolarity A is straightforward; we consider 
the case A = 3 in more detail below. 



III. TRANSITION STRENGTHS IN THE SPHERICAL LIMIT 

Close to the spherical limit, the deformation parameters of the intrinsic wave function are 
small and we can expand \<p(q)) to second order in q 2fl . The wave function is then projected 
on good angular momentum using the projection operator 

2 7+1/" 

P J mk = -g±- J dQV J MK (n)Rn (4) 

and the transformation properties of the Z operators. The ground state |0 + ) is obtained by 
projecting with Pq . It is given up to second order in q 2fl by 

|0 + ) = Ao {|0(O)> + ql {[Z ® Z\l + \z'^ |0(O)> + . . .} (5) 

Here we have introduced the notation q\ = ^^^^(-l) 2 ^ and 

[Z ® Z\ J M = J2(^'\JM)Z 2 ^Z 2 ^ (6) 
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Only the first term in Eq. (1), zeroth order in q 2jfl , will be required in the derivations below. 
The projection on J = 2 with the operator P^ IM gives the excited |2 + M) state as 

|2+M) = Af 2M {(-1) A Vm^2m|0(O)> + 0{q 2 2M )) (7) 

with a normalization factor H 2 m given by 

1 = IAWI 2 [ql M (m\Z^ M Z 2M \m) + 0{ql M )) . 

Since|0(O)) is a spherical wave function, the state Z 2 m\4>{0)) has angular momentum 2 and 
the mean value on the right hand side of the above equation is independent of M. It will 
be written as (||Z^Z 2 ||) which is a notation reminiscent of the reduced matrix elements of 
the Wigner-Eckart theorem. With this definition we finally obtain the expression for the 
normalized excited state wave function 

|2+M> = , Z l M - |0(O)) + 0(q 2M ) (8) 

(\\z+z 2 \\y/ 2 

The wave function |2 + M) is well defined in the q 2fl — > limit and is a linear combination of 
2-quasiparticle excitations of the spherical state. The expressions in Eqs (jSj) and can be 
now used in the defining formula for the B(E2) transition strength 

B(E2,0+ -> 2+) = JL |(2 + M|Q^|0 + )| 2 (9) 

AI/j, 

where Q\ is the standard electric multipole operator of rank A. Taking the expressions for 
the wave functions in the small deformation limit, the matrix element becomes 

(<p(0)\(z 2M ) + Q e 2fl \m) = ^M{\\ztQ e 2 \\) 

The final expression for the B(E2) is 

B( ^o^2 + ) |Sp ^ 5 AKij|g! (10) 

The generalization to arbitrary multipolarity A is 

B(£A ,0^A^) |Sph = (2 A +1 )^K^!. (u) 



In contrast to the rotational formula, Eq f lTOl) is nonzero in the spherical limit. This is a 
clear indication of the inadequacy of the rotational formula for the evaluation of transition 
strengths near sphericity. 



The quantities entering Eqs (TTO]) and (TTTj) can be calculated in linear response theory, but 
it is rather easy to calculate them using the intrinsic states of the HFB theory The only 
additional computational capability needed is the evaluation of matrix elements between 
different intrinsic states. In particular, we make use of the matrix element of quadrupole 
operator between deformed and spherical states given by 

(0fe,)IQLl0(O)) = -iq2A\\ZtQ e 2 \\) + 0(q 2 2u ). (12) 

To get the normalization in Eq. (7), we make use of the derivatives of the overlap function. 
The second derivative of the overlap between two intrinsic wave functions satisfies 

d 2 - - 

7 = T> — 7hT~ (0(?2i/)|0(gL'))|9 2 ?^o = (\\Z+Z 2 \\)5 UU/ . (13) 
0( l2uOq 2v , 

The second derivative can be approximated by a finite difference formula in the limit q 2v — > 



92^0 4gf„ 
Using this result and Eq ffl2l) we obtain the following result for the B(E2) in the spherical 
limit, 



5 K0WIQIJ0(O)}| 
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B(E2,0 + ->■ 2 + )| SDh = 5— lim 1 Hr m — — . (15) 

1 )lP 47r^oI(2-(0(g 2y )|0(-g 2y ))-(0(-g 2 ,)i0(g 2l/ )}) 1 1 

It is worth remarking that this derivation is valid for any value of v and therefore the axial 
case corresponding to v = can be used as well. This formula could be easily implemented 
in Wood-Saxon codes to obtain a quick estimate of the spherical transition strength. 

If this reasoning is applied to the octupole case, the |3~M) wave function is given by the 
expression 

13-M) = „ Z l M 10(0)) + 0(q 3M ). (16) 

This coincides with the negative parity projected wave function |\l/_(g3 M )) = A/"_(l — 
II)|0(g3 M )) up to order q^^. On the other hand, the |0 + ) wave function is given by the 
positive parity projected wave function |^/ + (g 3/i )) = A/+(l + n)|0(g 3A1 )) = |0(O)) + 0(q 2 ^). 
Taking into account these quantities in the general definition of Eq f TTTT) we arrive at the 
formula 

B(£3,0+ 3-), Sph « 7^|(^_(g3o)i^ i^ + (0)}| 2 (17) 
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Use of this formula of course requires that the ^30 in the negative parity wave function is small 
enough so that this wave function is well approximated by Eq ffT6]) . We used this formula 
recently in a global study of octupole correlations [9] to understand some discrepancies 
observed in the comparison with experimental data. 

We finish this section by mentioning that the previous methodology can also be used with 
scalar operators like the Hamiltonian. It is possible to obtain in this way formulas for the 
energies of J 7^ states in the spherical limit. This is briefly discussed in the appendix. 



IV. COMPARISON WITH EXACT PROJECTED TRANSITION STRENGTHS 



A. Validity of rotational formula 

In this section we compare the transition strengths computed with exact angular mo- 
mentum projection with the rotational formula and our spherical limit. The mean-field 
wave functions were calculated in the Hartree-Fock-Bogoliubov approximation assuming ax- 
ial symmetry and obtaining a range of deformations by including an external quadrupole 
field in the Hamiltonian. The range of deformations 02 spans the interval —0.3 to 0.4 in 



steps of 0.02 and a finer mesh with a step size of 0.01 is used in the -0.1 to 0.1 interval 



id- 



For those intrinsic wave functions the B(E2, 2 + — > 0) transition strength has been computed 
formula and exact angular momentum projection with 1 ( /?2 ) ) as the in- 



with the rotational 

trinsic states (see 11] for the relevant formulas). In Fig. [T]the ratio B(E2)hot / B(E2)proj 
is plotted as a function of 02 f° r a sample of nuclei spanning a wide range of masses. As 
expected, the ratio increases toward one as 02 becomes large. However, the limit is only 
reached in medium and heavy nuclei within our range of 2 values. For small values of 
02 the ratio is smaller than one and approaches zero as 2 — » 0. One can see that the 02 
value by itself does not provide a good indicator of the region of validity of the rotational 
formula. To get a more robust criterion, we go back to the basic assumption in deriving the 
rotational formula, that the intrinsic states have vanishing overlaps under finite rotations 
of the orientations. This requires a large angular momentum content of the intrinsic states. 
The mean square angular momentum of the intrinsic state (AJ 2 ) can be easily computed 
from the HFB wave function, so we may consider that quantity as a practical indicator. We 
note that overlap between rotated wave functions approaches a Gaussian of width 1/ (AT 2 ) 
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[3j. This result suggests that the validity of the rotational formula could be linked to spe- 
cific values of (AJ 2 ). To explore this possibility we have determined the value of (A J 2 ) 
for the intrinsic configuration that satisfies B(E2)phoj/B(E2)hqt ~ 3/4 (a value we have 
chosen to establish the limits of validity of the rotational formula) in each of the nuclei of 
our calculation. The values are shown as a histogram in Fig. [2J We see that the values are 
strongly peaked around (AJ 2 ) »s 10h 2 . This remarkable fact gives us an easily computed 
estimator of the validity of the rotational formula for the B(E2) transition strength for any 
nucleus in the Chart of Nuclides. 



B. Selected isotope and isotone chains 

The behavior of the spherical B(E2) transition strengths as a function of proton and 
neutron numbers is analyzed next. In Fig [3] the spherical transition strengths of Eq (Till]) 
are plotted as a function of neutron number for several isotopic chains. They have been 
computed using the exact angular momentum projected transition strengths for a deforma- 
tion of the intrinsic state of @2 = 0.005. The values of those spherical transition strengths 
are smaller than the typical values of well deformed nuclei that can reach a few hundreds 
of W.u. for heavy nuclei. The decrease with neutron number is rather weak except around 
magic neutron numbers where a marked peak is observed. This is probably a consequence of 
the lowering of the level density near magic numbers. Surprisingly, a peak at the non-magic 
number N = 40 is also seen. This behavior is not observed when the quantity is plotted 
as a function of proton number (see right panel). First the spherical transition strength 
increases with increasing Z values and a reduction at those values of Z corresponding to 
magic numbers is observed, specially at Z=82. The values of B(E2)s p h expressed in W.u. 
follow a trend with Z that is consistent with the expected linear behavior in Z based on 
the scaling of the mean value of proton's quadrupole moment (remember that W.u. scale 
like nuclear radius squared). A least square fit to the computed values for over two hundred 
nuclei yields the rule B(E2, 2 -> 0) Sph = 0.85 Z (W.u.). 
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Figure 1: The ratio B(E2)jiot / B(E2)p^oj 1S plotted as a function of the deformation parameter 
/?2 for a range of nuclei. The solid line connects calculated values. The dashed line is calculated 
from the interpolating formula, Eq. (|18|) and (JTUJ). 

C. An interpolating formula 



Even better than a criterion for the validity of Eq. (1) would be an interpolating formula 
that would also capture the transition region between spherical and strongly deformed nuclei. 
To this end we consider parameterizing the B(E2) by the function 

C 



£(£2,2+->0) Int 



exp[-(/V/^ 0) ) 2 ; 



■Pi 



(18) 
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Figure 2: Lowest (A. J 2 ) values of intrinsic wave functions that meet our criterion for using the 
rotational formula (see text). 

The parameter Cq is set to Co = (9e 2 ) / (80tt 2 )Z 2 Rq to recover the rotational formula at large 
deformation. The parameter is set to a value that reproduces the spherical limit, 

0®*= ^B{E2X -+0 + ), Sph . (19) 

The results obtained with Eq ([TBI) are plotted as dashed lines in Fig[TJ Remarkably, for most 
of the cases and for almost the whole range of 02 values both the exact and the approximate 
results are indistinguishable. It seems that our model can be used with confidence to compute 
B(E2) values provided that the parameter 0^ can be obtained. 

D. Computing the spherical limit 

An alternative formula for the evaluation of B(E2)s p h was obtained in Eq (fT5|) in terms 
of simple overlaps with the wave functions \4>{q2u))- To test its applicability we have per- 
formed calculations with our axially symmetric wave functions as a function of 2 an d some 
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Figure 3: The spherical-limit transition strengths of Eq (|10j) are displayed for several isotopic chains 
as a function of neutron number N on the left panel and as a function of Z in the right panel. The 
isotopic chains correspond to Z values between 12 and 94 in steps of 6 units. Strengths are given 
in Weisskopf units, 1 W.u. = 5.94 x 1(T 6 ,4 4 / 3 e 2 b 2 . 

representative results are given and compared to the exact results in Figure HI From the 
comparison we conclude that the formula is accurate enough for (3% values up to 0.05 for 
light nuclei and up to 0.01 for heavy ones and therefore can be used for a computationally 
inexpensive estimation of B(E2)s p h to be used in the model of Eq (Il8j) to compute the f3^ 
parameter as = (B(E2)s p h/Co) 1 ^ 2 ■ 

E. Octupole transitions 

Another interesting case to study is the one of the B(E3,3~ — > 0) transition strengths. 
They are associated to the octupole degree of freedom, parameterized in terms of the oc- 
tupole moments q^^. The rotational formula, valid in the strong quadrupole deformation 
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Figure 4: Accuracy of Eq. (|15p to calculate B{E2) values in the spherical limit. Plotted is the 
ratio of the B(E2) from Eq. (|15p to the value obtained by a full projection of the wave functions 
at /J2 = 0.005. The horizontal scale gives the /?2 values used in Eq. (I15p . The ratios are offset for 
clarity with the dotted lines indicating equal values. 

limit, reads in this case B(E3,3~ — > 0) = ^KQDl 2 - Contrary to the quadrupole deforma- 
tion case, there is no spontaneous parity symmetry breaking in most of the nuclei of the 
Nuclide chart with the exception of a few light Ra and Th isotopes and some rare earth 
nuclei like neutron poor Ba isotopes. Therefore the mean value of the octupole operator in 
the intrinsic state is zero. As a consequence, theories dealing with dynamical correlations 
are required in order to describe octupole correlations and the associated B{E3). In those 
theories the intrinsic octupole deformed state for the + is different from the one of the 3~. A 
typical example is that of parity projection with restricted variation of the intrinsic state 
that assigns the intrinsic states of the + and 3~ states to the ones producing the lowest par- 
ity projected energies E±(q 3 ) computed for axially symmetric octupole constrained intrinsic 
states with octupole deformation q 3Q . In this theory, the rotational formula restricted to axi- 
ally symmetric configurations becomes B(E3, 3~ — > + )|rot = ^\{^ -(Qsq^IQIoI^ +(q^))\ 2 
where now | ± (^30) ) are parity projected wave functions obtained from an intrinsic state 
with octupole deformation §30- In order to study the validity of this formula in the spherical 
limit, calculations as a function of the quadrupole moment should be carried out. The dif- 
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ficulty here is that there are two intrinsic states that potentially have different quadrupole 
deformations and therefore a study in terms of four variables (the quadrupole and oc- 
tupole moment of the positive and negative parity intrinsic states) should be carried out 
for a series of isotopes. Instead of this long calculation we have just taken the intrin- 
sic states for positive and negative parity from the results of |9| and computed the corre- 
sponding transition strengths with angular momentum projected wave functions. The ratio 
£(£3,3" -»■ + )| PRO j/ J B( J E3,3" -»■ + )| ROT is plotted in Fig [5] as a function of the 2 {+) 
deformation parameter of the positive parity intrinsic state. Values corresponding to nuclei 
where the negative parity quadrupole deformation parameter 02(~) differs significantly from 
(by ±0.1) have not been included in the plot. This includes nuclei with strong shape 
coexistence where the ground state is, for instance, prolate and the negative parity state 
is oblate. As a consequence of the mismatch in quadrupole deformations the overlap be- 
tween the wave functions is very small and the corresponding B(E3) are much smaller (and 
therefore more dependent on little details) than for intrinsic states with similar quadrupole 
deformation parameters. 

The first noteworthy observation is that the transition strengths computed with the 
projected angular momentum wave functions are always greater or equal the values obtained 
with the rotational formula. The results show that for 2 {+) values greater than 0.15 the 
rotational formula works reasonably well within a factor of 2. Around 02{+) = the 
ratio lies in between 3 and 8 in good agreement with the results of Eq (117]) that predict a 
factor 7 difference with the rotational formula in the spherical limit. The main conclusion 
is that for quadrupole deformations smaller than 2 ~ 0.15 the rotational formula should 
not be trusted and its use avoided in relating transition strengths to intrinsic octupole 
deformation parameters. A typical example illustrating the general trend is that of 208 Pb 
where the rotational formula predicts a B(E3) value of 7.1 W.u. whereas the transition 
strength with the angular momentum projected wave functions is 23.1 W.u. which is in 
much better agreement with the experimental data of 34 W.u.. Such enhancement of the 
B(E3) transition probabilities for near spherical configurations as compared to the rotational 



formula was already noticed in 
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13] for some spherical or near spherical nuclei. 
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Figure 5: The ratio B(E3,3 —> + )\prqj/B(E3, 3 — ^ + )|rot i s plotted as a function of the 
/%(+) quadrupole deformation parameter of the positive parity intrinsic state. 

V. CONCLUSIONS 

The validity of the rotational formula for multipole transition strengths is questioned for 
near spherical configurations. A general formula to compute those transitions in terms of 
intrinsic mean values and/or overlaps is derived by exploiting the simple structure of angular 
momentum projected wave functions in the spherical limit. An enhancement factor of 2A + 1 
for transitions of order A is obtained. Thorough numerical calculations of B(E2) and B(E3) 
transition strengths show the validity of the formulas obtained and establish criteria of 
validity for the rotational approximation. For quadrupole transitions, we proposed a simple 
model to compute the B(E2) and found that it is quite accurate over the entire range 
of deformation. The model contains two parameters that are fixed from the calculated 
transition strengths at the two limits, Eq. (2) and (9). only one parameter, ffl that 
unfortunately is nucleus dependent. We have also established a criteria to determine the 
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validity of the rotational formula that only requires the evaluation of a mean field quantity: 
the fluctuation (AJ 2 ) should be larger than ~ 10 for the rotational formula to be useful; 
it becomes quite accurate above (AJ 2 ) > 15. For octupole transition strengths B(E3), 
the quadrupole deformation parameter /3 2 of the ground state has to be larger than 0.15 
for the rotational formula to be valid and precautions are in order for those cases of shape 
coexistence where the quadrupole deformation parameters of positive an negative parity 
states differ considerably. For spherical configurations the B(E3) can be up to a factor of 8 
larger than the values provided by the rotational formula. 

A table is provided, as supplementary material, with the spherical B(E2) strengths and 
the (5^ parameters for 818 even-even nuclei computed with the Gogny D1S interaction. 
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Appendix A: Projected energies in the spherical limit 

The same arguments used in the previous section can be used to compute the energy of 
the |2 + M) as given by Eq (jSJ) in the spherical limit 

( w "liW +0( * ,J - 

Defining 

V = T> — ar^WI^I^DlB^-m = (\\Z^HZ 2 \\)6 vvl 
as using Eq ( TTBl the excitation energy can be written as 



h nn i 

Sph = 



E{2 + ) v 
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an expression that coincides with twice the zero point energy correction obtained in the 
Generator Coordinate Method (GCM) for the quadrupole coordinate in the harmonic limit 
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of the Gaussian Overlap Approximation (GOA) (See Eq (10.136) of [3]). The energy of the 
2 + state in the spherical limit is not given in calculations with angular momentum projection 
as its evaluation involves the ratio of two very small quantities which are difficult to compute 



with the required accuracy 
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